Phylogenetic relations and range history of jerboas of the Allactaginae subfamily (Dipodidae, Rodentia).

Five-toed jerboas of the subfamily Allactaginae comprise several complex taxa occurring over a wide distribution range covering a large part of the Eurasian arid belt. In this study, we employed current methods of molecular phylogenetics based on 15 nuclear genes and the mitochondrial gene cytb to revise relations and systematics within Allactaginae. We also applied species distribution modelling projected on paleo-environmental data to reconstruct the geographic patterns of speciation in Allactaginae. We elucidated the intergeneric relationships within this subfamily and clarified interspecies relations within the genus Scarturus. Moreover, our results demonstrate the species status of S. caprimulga; outline the currently understudied diversity within Orientallactaga, Allactaga, and Pygeretmus; and improve the divergence estimates of these taxa. Based on our results from modelling of geographic range fragmentation in allactagines, we suggest the dating and location of speciation events and present hypotheses regarding general habitat niche conservatism in small mammals.

An Early Oligocene age for the oldest known monkeys and rodents of South America.

The Santa Rosa fossil locality in eastern Perú produced the first Paleogene vertebrate fauna from the Amazon Basin, including the oldest known monkeys from South America. This diverse paleofauna was originally assigned an Eocene age based largely on the stage of evolution of the site's caviomorph rodents and marsupials. Here, we present detrital zircon dates that indicate that the maximum composite age of Santa Rosa is 29.6 ± 0.08 Ma (Lower Oligocene), although several zircons from Santa Rosa date to the Upper Oligocene. The first appearance datum for Caviomorpha in South America is purported to be the CTA-27 site in the Contamana region of Perú, which is hypothesized to be ∼41 Ma (Middle Eocene) in age. However, the presence of the same caviomorph species and/or genera at both CTA-27 and at Santa Rosa is now difficult to reconcile with a >11-My age difference. To further test the Middle Eocene age estimate for CTA-27, we ran multiple Bayesian tip-dating analyses of Caviomorpha, treating the ages of all Paleogene species from Perú as unknown. These analyses produced mean age estimates for Santa Rosa that closely approximate the maximum 29.6 ± 0.08 Ma composite date provided by detrital zircons, but predict that CTA-27 is much younger than currently thought (∼30 Ma). We conclude that the ∼41 Ma age proposed for CTA-27 is incorrect, and that there are currently no compelling Eocene records of either rodents or primates in the known fossil record of South America.

Prevalence of Hantaviruses Harbored by Murid Rodents in Northwestern Ukraine and Discovery of a Novel Puumala Virus Strain.

In Europe, two species of hantaviruses, Puumala orthohantavirus (PUUV) and Dobrava orthohantavirus (DOBV), cause hemorrhagic fever with renal syndrome in humans. The rodent reservoirs for these viruses are common throughout Ukraine, and hence, the goal of this study was to identify the species and strains of hantaviruses circulating in this region. We conducted surveillance of small rodent populations in a rural region in northwestern Ukraine approximately 30 km from Poland. From the 424 small mammals captured, we identified nine species, of which the most abundant were Myodes glareolus, the bank vole (45%); Apodemus flavicollis, the yellow-necked mouse (29%); and Apodemus agrarius, the striped field mouse (14.6%) Using an indirect immunofluorescence assay, 15.7%, 20.5%, and 33.9% of the sera from M. glareolus, A. glareolus, and A. flavicollis were positive for hantaviral antibodies, respectively. Additionally, we detected antibodies to the hantaviral antigen in one Microtus arvalis, one Mus musculus, and one Sorex minutus. We screened the lung tissue for hantaviral RNA using next-generation sequencing and identified PUUV sequences in 25 small mammals, including 23 M. glareolus, 1 M. musculus, and 1 A. flavicollis, but we were unable to detect DOBV sequences in any of our A. agrarius specimens. The percent identity matrix and Bayesian phylogenetic analyses of the S-segment of PUUV from 14 M. glareolus lungs suggest the highest similarity (92-95% nucleotide or 99-100% amino acid) with the Latvian lineage. This new genetic information will contribute to future molecular surveillance of human cases in Ukraine.

Assessment of crop damage by rodent pests from experimental barley crop fields in Farta District, South Gondar, Ethiopia.

This study was conducted in Farta district, south Gondar from 2019 to 2020 cropping years to identify rodent pest species and estimate damage caused on barley crops. Four independent barley crop fields (40 x 40 m each) were sampled randomly to estimate the loss. Two were located near Alemsaga Priority State Forest and the other two were away from the forest. Four (2 x 2 m) rodent exclusion plots were established at 10 m interval as control units in each selected experimental barley fields using fine wire mesh. Rodent pest species were collected using both Sherman and snap traps throughout the different crop growing stages. The damaged and undamaged barley tillers by pest rodents were counted on five 1 x 1 m randomly sampled quadrats for each selected experimental fields. Variations on pest rodent population between cropping years and sites were analyzed using Chi square test. The mean crop damages between cropping years and experimental field sites were analyzed using two way ANOVA. Arvicanthis abyssinicus, Mastomys natalensis, Arvicanthis dembeensis, Mus musculus, Lophuromys simensis, Tachyoryctes splendens and Hystrix cristata were identified as pest rodents in the study area. A total of 968 individual rodents (427 in 2019 and 541 in 2020) were trapped during the study period. There was a statistical variation (χ2 = 13.42, df = 1 and P<0.05) between trapped individuals of the two successive years. The crop fields near the forest were more vulnerable than away from the forest during both cropping years. Statistical variations was observed on mean crop losses between cropping years and experimental barley crop sites. The highest crop damage was seen at maturity stage and the lowest during sowing in all experimental plots and cropping years. The percentage of barley yield loss due to rodent pests was 21.7 kg ha-1. The monetary value of this yield loss was equivalent to 4875 Birr (121.9 US$ h-1). Alemsaga Forest as shelter and conservation strategies like free of farmland from livestock and terracing for soil conservation have great role for the high rodent pest populations in the study area. Field sanitation, trapping and using restricted rodenticides like zinc phosphide are the possible recommendation to local farmers against rodent pests.

Small mammal glucocorticoid concentrations vary with forest fragment size, trap type, and mammal taxa in the Interior Atlantic Forest.

Species that live in degraded habitats often show signs of physiological stress. Glucocorticoid hormones (e.g., corticosterone and cortisol) are often assessed as a proxy of the extent of physiological stress an animal has experienced. Our goal was to quantify glucocorticoids in free-ranging small mammals in fragments of Interior Atlantic Forest. We extracted glucocorticoids from fur samples of 106 small mammals (rodent genera Akodon and Oligoryzomys, and marsupial genera Gracilinanus and Marmosa) from six forest fragments (2-1200 ha) in the Reserva Natural Tapytá, Caazapá Department, Paraguay. To our knowledge, this is the first publication of corticosterone and cortisol levels for three of the four sampled genera (Akodon, Oligoryzomys, and Marmosa) in this forest system. We discovered three notable results. First, as predicted, glucocorticoid levels were higher in individuals living withing small forest fragments. Second, animals captured live using restraint trapping methods (Sherman traps) had higher glucocorticoid levels than those animals captured using kill traps (Victor traps), suggesting that hair glucocorticoid measures can reflect acute stress levels in addition to long-term glucocorticoid incorporation. These acute levels are likely due to urinary steroids diffusing into the hair shaft. This finding raises a concern about the use of certain trapping techniques in association with fur hormone analysis. Finally, as expected, we also detected genus-specific differences in glucocorticoid levels, as well as cortisol/corticosterone ratios.

Rickettsia spp. in rodent-attached ticks in Estonia and first evidence of spotted fever group Rickettsia species Candidatus Rickettsia uralica in Europe.

BACKGROUND: Rickettsia spp. are human pathogens that cause a number of diseases and are transmitted by arthropods, such as ixodid ticks. Estonia is one of few regions where the distribution area of two medically important tick species, Ixodes persulcatus and I. ricinus, overlaps. The nidicolous rodent-associated Ixodes trianguliceps has also recently been shown to be present in Estonia. Although no data are available on human disease(s) caused by tick-borne Rickettsia spp. in Estonia, the presence of three Rickettsia species in non-nidicolous ticks has been previously reported. The aim of this study was to detect, identify and partially characterize Rickettsia species in nidicolous and non-nidicolous ticks attached to rodents in Estonia. RESULTS: Larvae and nymphs of I. ricinus (n = 1004), I. persulcatus (n = 75) and I. trianguliceps (n = 117), all removed from rodents and shrews caught in different parts of Estonia, were studied for the presence of Rickettsia spp. by nested PCR. Ticks were collected from 314 small animals of five species [Myodes glareolus (bank voles), Apodemus flavicollis (yellow necked mice), A. agrarius (striped field mice), Microtus subterranius (pine voles) and Sorex araneus (common shrews)]. Rickettsial DNA was detected in 8.7% (103/1186) of the studied ticks. In addition to identifying R. helvetica, which had been previously found in questing ticks, we report here the first time that the recently described I. trianguliceps-associated Candidatus Rickettsia uralica has been identified west of the Ural Mountains.

Molecular and morphological evidence for a new species of the genus Typhlomys (Rodentia: Platacanthomyidae).

In this study, we reassessed the taxonomic position of Typhlomys (Rodentia: Platacanthomyidae) from Huangshan, Anhui, China, based on morphological and molecular evidence. Results suggested that Typhlomys is comprised of up to six species, including four currently recognized species ( Typhlomys cinereus, T. chapensis, T. daloushanensis, and T. nanus), one unconfirmed candidate species, and one new species ( Typhlomys huangshanensis sp. nov.). Morphological analyses further supported the designation of the Huangshan specimens found at mid-elevations in the southern Huangshan Mountains (600 m to 1 200 m a.s.l.) as a new species.

Adaptive evolution of ß-globin gene in subterranean in South America octodontid rodents.

The convergent evolution of subterranean rodents is an excellent model to study how natural selection operates and the genetic bases of these adaptations, but the study on the different taxa has been very uneven and still insufficient. In the octodontoid caviomorph rodent superfamily there are two independent lineages where they have recently evolved into totally underground lifestyles: the genera Ctenomys (tuco-tucos) and Spalacopus (coruro). The underground habitat is characterized by an hypoxic and hypercapnic atmosphere, thus gas exchange is one of the most important challenges for these animals. The invasion of the underground niche could have modified the selective regimes of proteins involved in the respiration and transport of O2 of these rodents, positively selecting mutations of higher affinity for O2. Here we examine the sequence variation in the beta globin gene in these two lineages, within a robust phylogenetic context. Using different approaches (classical and Bayesian maximum likelihood (PAML/Datamonkey) and alternatives methods (TreeSAAP)) we found at least three sites with evidence of positive selection in underground lineages, especially the basal branch that leads to the Octodontidae family and the branch that leads to the coruro, suggesting some adaptive changes to the underground life. We also found a convergence with another underground rodent, which cannot be identified by the above methods.

Orthohantaviruses infections in humans and rodents in Baoji, China.

In recent years, hemorrhagic fever with renal syndrome (HFRS) incidence has been becoming a severe public health problem again due to its significant increase in Shaanxi Province, China. Baoji, located in the Guanzhong Plain in the central part of Shaanxi Province, has been severely affected by HFRS since its first emergence in 1955. To better understand the epidemiology of orthohantaviruses infection in humans and the causative agents carried by the rodents, the long-term incidence patterns were analyzed and a molecular epidemiological investigation of orthohantaviruses infection in humans and rodents was performed. During 1984-2019, 13,042 HFRS cases were registered in Baoji, including 275 death cases. Except the first high prevalence of HFRS in 1988-1993, another two epidemic peaks were observed in 1998-2003 and 2012, respectively, although vaccination project was started since 1996. During the same period, HFRS cases in Baoji mainly were recorded in winter suggesting they may be caused by Hantaan orthohantavirus (HTNV), while a small peak of HFRS was also found in summer with unknown reason. Nucleotide identity and phylogenetic analyses demonstrated that a novel clade of HTNV sequences recovered from HFRS cases were closely related to those from rodents, including species close contact with humans, suggesting a direct viral transmission from rodents to humans and the important role in the HTNV transmission the nontraditional rodent hosts may play. Moreover, two distant related Dabieshan orthohantavirus (DBSV) lineages were also identified in Niviventer niviventer in this area demonstrating its considerable genetic diversity. Our data indicated that continual spillover of HTNV from rodents to humans, contributing to the high prevalence of HFRS in humans in Baoji.

Scalable classification of organisms into a taxonomy using hierarchical supervised learners.

Accurately identifying organisms based on their partially available genetic material is an important task to explore the phylogenetic diversity in an environment. Specific fragments in the DNA sequence of a living organism have been defined as DNA barcodes and can be used as markers to identify species efficiently and effectively. The existing DNA barcode-based classification approaches suffer from three major issues: (i) most of them assume that the classification is done within a given taxonomic class and/or input sequences are pre-aligned, (ii) highly performing classifiers, such as SVM, cannot scale to large taxonomies due to high memory requirements, (iii) mutations and noise in input DNA sequences greatly reduce the taxonomic classification score. In order to address these issues, we propose a multi-level hierarchical classifier framework to automatically assign taxonomy labels to DNA sequences. We utilize an alignment-free approach called spectrum kernel method for feature extraction. We build a proof-of-concept hierarchical classifier with two levels, and evaluated it on real DNA sequence data from barcode of life data systems. We demonstrate that the proposed framework provides higher f1-score than regular classifiers. Besides, hierarchical framework scales better to large datasets enabling researchers to employ classifiers with high classification performance and high memory requirement on large datasets. Furthermore, we show that the proposed framework is more robust to mutations and noise in sequence data than the non-hierarchical classifiers.

Karyotypic diversity within the genus Makalata (Echimyidae: Echimyinae) of Brazilian Amazon: Chromosomal evidence for multiple species.

The genus Makalata is a taxonomically complex group of rodents on which few cytogenetic studies have been performed. Most of the published karyotypes were described based only on conventional chromosome staining. Here, we studied the karyotypes of Makalata from two Brazilian Amazonian states, Amapá and Pará, by Giemsa-staining, G- and C-banding, AgNO3-staining and FISH with 18S rDNA and telomeric sequences probes. We observed 2n = 66/FN = 124 in the Pará state population in Makalata sp; and 2n = 72/FN = 128 in the Amapá state population in M. didelphoides. Multiple chromosome rearrangements may have given rise to these karyotypes, which differ significantly from each other and from those reported in the literature. The chromosomal differences among the described Makalata karyotypes can act as a barrier to gene flow; since they are also associated with geographic barriers (e.g., rivers) and numerous molecular differences, they could be seen as evidence for reproductive isolation of populations from genus Makalata. Our data suggest that the genus is chromosomally diverse and the karyotypes may belong to different species. These karyotypes may prove useful as taxonomic markers for these rodents.

A Retrospective Survey of Rodent-borne Viruses in Rural Populations of Brazilian Amazon.

INTRODUCTION: The Amazon tropical rainforest has the most dense and diverse ecosystem worldwide. A few studies have addressed rodent-borne diseases as potential hazards to humans in this region. METHODS: A retrospective survey was conducted using enzyme-linked immunosorbent assay for detecting mammarenavirus and orthohantavirus antibodies in 206 samples collected from rural settlers of the Brazilian Western Amazonian region. RESULTS: Six (2.91%) individuals in the age group of 16 to 36 years were found to possess antibodies against mammarenavirus. CONCLUSION: Evidence of previous exposure to mammarenavirus in the rural population points to its silent circulation in this region.

Evolutionary analysis of chemokine CXCL16 and its receptor CXCR6 in murine rodents.

The chemokine CXCL16 and its receptor CXCR6 are implicated in various physiological and pathological processes in cooperative and/or stand-alone fashions. Despite the significance of rodent animal models in elucidating the function and clinical relevance of the chemokine and its receptor, the evolutionary characterization of these molecules remains deficient for this taxon to a certain extent. In this study, we implemented a comparison of synonymous and nonsynonymous variation rates in combination with the maximum likelihood (ML) analysis and Tajima's test to evaluate the interspecific and intraspecific evolutions of CXCL16 and CXCR6 in murine rodents. Our results indicate that adaptive selection has frequently contributed to genetic diversity of both CXCL16 and CXCR6 in the murine lineage that is asynchronous with a relative dependence between these genes. This signature is radically different from the lineage-specific and concordant adaptive diversity of the primate homologues of these genes, which was reported in a previous study. The diversity identified in the present study shed further light on molecular strategies against the challenges towards CXCL16 and CXCR6.

A parapithecid stem anthropoid of African origin in the Paleogene of South America.

Phylogenetic evidence suggests that platyrrhine (or New World) monkeys and caviomorph rodents of the Western Hemisphere derive from source groups from the Eocene of Afro-Arabia, a landmass that was ~1500 to 2000 kilometers east of South America during the late Paleogene. Here, we report evidence for a third mammalian lineage of African origin in the Paleogene of South America-a newly discovered genus and species of parapithecid anthropoid primate from Santa Rosa in Amazonian Perú. Bayesian clock-based phylogenetic analysis nests this genus (Ucayalipithecus) deep within the otherwise Afro-Arabian clade Parapithecoidea and indicates that transatlantic rafting of the lineage leading to Ucayalipithecus likely took place between ~35 and ~32 million years ago, a dispersal window that includes the major worldwide drop in sea level that occurred near the Eocene-Oligocene boundary.

Seroprevalence of leptospiral antibodies in rodents from riverside communities of Santa Fe, Argentina.

BACKGROUND: Leptospirosis is a zoonotic disease that can be transmitted by contact with the urine of infected mammals. Rodents play a mayor role in the transmission of leptospires to humans. The province of Santa Fe reports the greatest number of cases in Argentina. Yet, in this region, there are still knowledge gaps regarding the diversity of rodent species that may be hosts of pathogenic leptospires. The aims of this study were to evaluate the presence of leptospiral antibodies in rodents from three riverside communities of Santa Fe, and to identify factors associated with leptospiral infection. METHODOLOGY/PRINCIPAL FINDINGS: Each community was divided into three environmental settings based on the level of human disturbance, and sampled during two springs (Sep-Oct 2014 and 2015) and one autumn (Mar-Apr 2015). Serum samples of captured sigmodontine and murine rodents were tested for leptospiral antibodies by enzyme-linked immunosorbent assay (ELISA), and microagglutination test (MAT) was used to assess the infecting serovar in seropositive individuals. Factors influencing seropositivity were analyzed using logistic regression models. We caught 119 rodents, of which 101 serums were suitable for analysis. Most frequently trapped species were Scapteromys aquaticus, Akodon azarae and Oligoryzomys spp., with seroprevalences of 41.3%, 42.9% and 55% respectively. Seropositivity was higher in individuals with an average body condition score and in those that were sexually mature, but in the latter the differences were marginally significant. CONCLUSIONS/SIGNIFICANCE: Our results suggest that native rodents may be playing a role in the environmental circulation of pathogenic leptospires and provide relevant information for public health policies in the area.

Chromosomal Signatures Corroborate the Phylogenetic Relationships within Akodontini (Rodentia, Sigmodontinae).

Comparative chromosome-painting analysis among highly rearranged karyotypes of Sigmodontinae rodents (Rodentia, Cricetidae) detects conserved syntenic blocks, which are proposed as chromosomal signatures and can be used as phylogenetic markers. In the Akodontini tribe, the molecular topology (Cytb and/or IRBP) shows five low-supported clades (divisions: "Akodon", "Bibimys", "Blarinomys", "Oxymycterus", and "Scapteromys") within two high-supported major clades (clade A: "Akodon", "Bibimys", and "Oxymycterus"; clade B: "Blarinomys" and "Scapteromys"). Here, we examine the chromosomal signatures of the Akodontini tribe by using Hylaeamys megacephalus (HME) probes to study the karyotypes of Oxymycterus amazonicus (2n = 54, FN = 64) and Blarinomys breviceps (2n = 28, FN = 50), and compare these data with those from other taxa investigated using the same set of probes. We strategically employ the chromosomal signatures to elucidate phylogenetic relationships among the Akodontini. When we follow the evolution of chromosomal signature states, we find that the cytogenetic data corroborate the current molecular relationships in clade A nodes. We discuss the distinct events that caused karyotypic variability in the Oxymycterus and Blarinomys genera. In addition, we propose that Blarinomys may constitute a species complex, and that the taxonomy should be revised to better delimit the geographical boundaries and their taxonomic status.

Review of Rodent Euthanasia Methods.

The optimal choice of euthanasia method for laboratory rodents depends on a number of factors, including the scientific goals of the study, the need to minimize animal pain and/or distress, applicable guidelines and laws, the training and proficiency of personnel, and the safety and emotional needs of the personnel performing the euthanasia. This manuscript aims to provide guidance to researchers so they may select the method of euthanasia that results in minimal experimental confounds, such as the creation of artifact and alteration of tissues and analytes. Specific situations addressed include euthanasia of large numbers of rodents and euthanasia of neonates. Recent literature supports the notion of significant strain-dependent differences in response to euthanasia methods such as CO2 inhalation. To assist researchers in selecting a strain-appropriate method of euthanasia, the authors present a summary of methodologies for assessing the effectiveness of euthanasia techniques, including elements and parameters for a scoring rubric to assess them.

Burrow systems evince non-solitary geomyid rodents from the Paleogene of southern Mexico.

We describe a new complex burrow system produced by geomyids in southern Mexico. Yaviichnus inyooensis igen. isp. nov. is composed of main large chambers near the top of the paleosol, from which shafts showing different morphologies and orientations radiate, some of them ending in or connected to small deeper chambers. Gregorymys spp. is proposed as the producer based on its fossorial habits, abundance in the outcrops, presence of remains inside the burrows, and paired grooves in the walls, which are compatible with the traces of geomyid incisors. The complexity of these burrows attests to an extended underground life that would have been triggered by semiarid to arid conditions. Morphological complexity also suggests that the burrows were excavated and inhabited by more than one individual, indicating that Oligocene Gregorymys of southern Mexico would be a unique gregarious geomyid.

Rodents as intermediate hosts of cestode parasites of mammalian carnivores and birds of prey in Poland, with the first data on the life-cycle of Mesocestoides melesi.

BACKGROUND: Rodents constitute an important part of the diet of many carnivore species. This predator-prey food chain is exploited by helminth parasites, such as cestodes, whose larval stages develop in rodents and then mature to the adult stage in predators. The main aim of our study was to use molecular techniques for identification of cestode species recovered from both intermediate and definitive hosts, with a particular focus on the genus Mesocestoides. METHODS: Larval cestodes were obtained during our long-term studies on rodent helminth communities in the Mazury Lake District in the north-east Poland in 2000-2018. Cestode larvae/cysts were collected from body cavities or internal organs (e.g. liver) during autopsies. Adult tapeworms were derived from nine red foxes, three Eurasian badgers and one Eurasian lynx. PCR amplification, sequencing and phylogenetic analyses were conducted employing three genetic markers: 18S rDNA, mitochondrial (mt) 12S rDNA and the mt cytochrome c oxydase subunit 1 (cox1) gene fragment. RESULTS: Altogether 19 Mesocestoides samples were analyzed, including 13 adult tapeworms from definitive hosts and six larval samples from 4 bank voles and 2 yellow-necked mice. Phylogenetic analyses revealed three well-supported trees of similar topology. In each case the Mesocestoides samples formed two separate clades. All isolates from foxes, the lynx isolate and two isolates from rodents grouped with Mesocestoides litteratus. Four isolates from rodents and all three isolates from Eurasian badgers were resolved in a separate clade, most similar to North American M. vogae (syn. M. corti). Examination of fixed, stained adult specimens from Eurasian badgers revealed consistency with the morphology of Mesocestoides melesi. Therefore, this clade is likely to represent M. melesi, a species first described in 1985 from the Eurasian badger Meles meles. Molecular analysis allowed also the identification of Taenia crassiceps, Hydatigera kamiyai and Cladotaenia globifera among larvae derived from rodents. CONCLUSIONS: Molecular and phylogenetic analyses support the recognition of M. melesi as a valid species. Our data represent the first record of the larvae of this species in rodents. This is the first report on the occurrence of H. kamiyai in rodents from Poland.

Evolutionary history of two cryptic species of northern African jerboas.

BACKGROUND: Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence. RESULTS: The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization. CONCLUSIONS: Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.